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ABSTRACT 

Aims. In this paper we present the first results of a detailed spectroscopic and photometric analysis of the V = 1 1.7™ eclipsing binary 
ASAS J052821+0338.5. 

Methods. With the FIES spectrograph at the Nordic Optical Telescope we have obtained a series of high-resolution spectra 
(R » 47000) covering the entire orbit of the system. In addition we obtained simultaneous broadband photometry from three small 
aperture telescopes. From these spectroscopic and photometric data we have derived the system's orbital parameters and determined 
the fundamental stellar parameters of the two components. 

Results. Our results indicate that ASAS J05282 1+0338.5 is a K1/K3 pre-main-sequence eclipsing binary, with component masses 
of 1.38Mq and 1.33 Mq and a period of 3.87 days, located at a distance of 280 ± 30 pc. The kinematics, physical location and the 
evolutionary status of the two stars suggest that ASAS J05282I+0338.5 is a member of the ~ 1 1 Myr old Orion OB la subassociation. 
The systems also exhibits smooth ~ 0.15'" out-of-eclipse variations that are similar to those found in RS CVn binaries. Furthermore 
the parameters we derived are consistent with the 10-13 Myr isochrones of the popular Baraffe stellar evolutionary models. 

Key words, stars: pre-main sequence - stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: fundamental parameters - 
stars: individual: ASAS J052821+0338.5 



1. Introduction 

The advantageous geometry of detached double-line eclipsing 
binary systems makes it possible to measure stellar radii and 
masses to very high accuracy (1-2%, Andersen 1991 ). Such ac- 
curate parameters are of crucial importance in providing ob- 
servational constraints on stellar evolutionary models. However, 
very few of the known double-lined eclipsing binaries are stars 
at early stages of evolution, leaving the models in this domain 
with only few empirical points of reference. 

The list of known double-lined eclipsing binaries with com- 
ponents that are in the pre-main-sequence (PMS) phase of evo- 
lution is not long. Up to now, only five low-mass double-lined 
eclipsing binaries in which both components are PMS objects 
have been reported, RXJ 0529.4+0041 A (Covino et al. iMjQl 
Covino et al. ,2004) . VI 174 Ori (Stassun et al. 2004), the brown 
d warf eclipsing binary 2MASS J05352184-0546085 (St assun et 
aL 120061 12OO7 I) and, very rec ently, J W 380 (Irwin et al. ! 2007bb 
and Par 1802 (Cargile et al. |20D71l. In addition, EK C ep and 
TY C iA (Popper 1980; Popper 1987; Andersen jlMTl Casey 
11998) are known to harbour one PMS component and one main- 
sequence component. 

In this paper, we announce the discovery that 
ASAS J052821 +0338.5 (also known as GSC2.2 N3003 11274, 
2MASS 05282082+0338327 and USNO-B 1.0 0936-0073790) 
is a PMS eclipsing binary system. In addition, the sys- 
tem could be associated with the ROSAT X-ray source 
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IRXS J052820.4+033823, located only 11" away, within 
ROSAT's positional error margin of 14". Using new high- 
cadence photometric and high-resolution spectroscopic 
observations we have determined the system's fundamental 
properties and confirm the PMS nature of both components. 



2. Observations 

2.1. Photometry 

ASAS J052821+0338.5 was first iden tified in the All-Sky 
Automated Survey (ASAS, Pojmanski "20021 as a detached 
eclipsing binary system with a maximum, out-of-eclipse V-band 
magnitude of V = 11 .7™ and a period of P = 1 .9364 days. In ad- 
dition, the target is listed in several large photometric databases 
consolidated in the Naval Observatory Merged Astronomical 
Dataset (NOMAD- 1.0, Zacharias et al. |2004ji, which provide 
optical and infra-red magnitudes of B = 12.68", J = 9.636™, 
H = 9.098" and K = 8.961". 

In the ASAS survey, 171 V-band measurements of the ob- 
ject were obtained during the period February 2001 -April 2006 
with a median sampling rate of 0.25 d"' . The resulting lightcurve 
exhibits periodic eclipses with a depth of ~ 0.3" and scatter 
in the out-of-eclipse phases with an rms amplitude of 0.15", 
likely due to an evolving spot configuration. We supplemented 
the long-term photometry from ASAS with additional precise, 
high-cadence observations obtained with several small aperture 
telescopes. 
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Fig.l. All differential photometi-y of ASAS J052821 +0338.5 
collected by the NMSU-lm, JGT-lm and 3MT, phase-folded 
with the orbital period of 3.8729 days and epoch. To - 
2454110.8500. The error bars on the data are smaller than the 
plotted symbols, and estimates of the errors are given in the text. 



2.1.1. High cadence, high precision photometry 

The New Mexico State University 1-m telescope (NMSU-lm) 
provided the majority of the follow-up photometry. A total of 
8805 V-band measurements of the target were acquired over 
31 nights between January 10-April 17, 2007. The data were 
processed with standard data reduction procedures including 2D 
bias and overscan subtraction, flat-fielding, and aperture pho- 
tometry. To derive the differential photometry lightcurve, the 
magnitudes of three nearby, bright stars were average combined 
into a non-variable reference object which was subtracted from 
the instrumental magnitude of the target star. On February 4, 
5, and 7, three additional nights of V-band photometry were 
obtained with the 1-m James Gregory Telescope (JGT-lm, 
St. Andrews, UK). These data can be directly compared with the 
NMSU-lm data as the same basic reduction procedures were 
applied to the images and the same three reference stars were 
used to derive the differential photometry. The V-band differen- 
tial photometry is normalized so that the magnitude of the phase 
between 0.94-0.96 (before the primary eclipse) has a median 
value of zero. 

Time-series photometry in a redder waveband was obtained 
with the 0.5m 3-Mirror Telescope (3MT; Cambridge, UK) from 
January-March, 2007. The filter (RG630) plus the CCD efli- 
ciency curve for the detector selects a broadband wavelength 
region which is between the Johnson R and Cousins /-band fil- 
ters. These data were reduced in a standard way with a version 
of the Cambridge Astronomical Survey Unit (CASU) data re- 
duction and photometry pipeline (Irwin & Lewis 120011 Irwin et 
al. 2007a) modified for the 3MT. These data are normalized to 
match the V-band photometry, so that the phases before the pri- 
mary eclipse have a median value of zero. 

The final differential photometry lightcurve of all telescopes 
folded with the ephemeris derived below is shown in Fig. [T] 
In Fig. |2] we show the original ASAS photometry (top panel) 
as well as our small-telescope data, divided into panels with 
photometry obtained within consecutive 2-4 week time inter- 
vals, the typical timescale for a stable starspot configuration. 
The ASAS photometry has a precision of ~ 0.05™. The small- 
telescope V-band photometry (crosses) has a sampling of < 
1 minute, and has an uncertainty ranging from ~ 0.003 to 
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Fig. 2. This figure illustrates the evolution of out-of-eclipse 
variabiHty on ASAS J052821+0338.5 during the four months 
(January- April, 2007) following the spectroscopic observations. 
The original ASAS photometry is shown in the top panel. 
Subsequent panels show the V-band lightcurve obtained with the 
NMSU-lm and JGT-lm telescopes (black crosses) as well as the 
red (RG630) lightcurve obtained with the 3MT (grey diamonds). 
The data are broken into 2-A week intervals. The time span (in 
HJD - 2450000) covered by the data in each panel is given in the 
top right hand corner The solid line is a simple one-spot model, 
plotted as a reference. 



0.010™, mainly due to variable systematic errors such as chang- 
ing sky transparency. We opted to exclude 17 measurements 
which had unusually large photometric errors of > 0.05 mag. 
The red wavelength data from the 3MT (diamonds) has a preci- 
sion of 0.006-0.008™ and has a slightly higher scatter and lower 
sampling due to the smaller aperture telescope used. 

The lightcurve shows deep primary and secondary eclipses 
(~ 0.37™), a period (P ^ 3.87 days) twice what was determined 
from the ASAS photometry alone, and large-amplitude out-of- 
eclipse variability (~ 0.15™). The out-of-eclipse variations are 
present for the duration of our four month observing campaign, 
but they are not completely stable. There are small variations in 
the lightcurve and in the depth of the secondary eclipse (approx- 
imately ±0.01™), probably caused by an evolving spot config- 
uration (see also the discussion in Sect. [4|. The observed out- 
of-eclipse variability is qualitatively consistent with a large cool 
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Fig. 3. This spectrum of ASAS J052821+0338.5 was obtained at phase 0.64, and shows clearly the double-lined nature of the 
spectrum, the emission in Ha as well as the strong absorption in the Li i 6708 line. 



spotted region on the back side of the primary star causing the 
attenuation of the secondary eclipse. Also, the secondary eclipse 
in the ASAS data appears to be only half the depth that it exhibits 
in the NMSU-lm data. This is a further suggestion that the to- 
tal depth of the secondary eclipse is being affected by variable 
star spots. 

2.1.2. Ephemeris 

Partial primary and secondary eclipses which were detected in 
the NMSU-lm data were used in combination with the ASAS 
photometry to derive an ephemeris for the system. We defined 
the epoch of the system, Tq, to be the midpoint of the most com- 
plete primary eclipse. We used the NMSU-lm data obtained on 
HJD=24541 10 which contains the only highly time sampled pri- 
mary eclipse showing the eclipse minimum to determine Tq. We 
modelled the shape of this eclipse with a 4th order polynomial 
and performed a minimization with respect to the model to 
derive the best fitting coefficents. We then set Tq to be the mid- 
point of the model curve. Next, we fixed To and solved for the 
period of the system using the nearly complete primary eclipse 
obtained on HJD=2454110 and the only other partial primary 
eclipse (HJD =24541 16) in the NMSU-lm data in combination 
with the complete ASAS lightcurve. We tested a grid of possible 
periods, phase-folding these data on each trial period, and found 
the period which minimized the of these data with respect to 
the model. The resulting ephemeris at the midpoint of the eclipse 
is given by 

T(V^in) = HJD 24541 10.8500 ± 0.002 + 3.8729£ + 0.0004 

where E corresponds to the cycle number The uncertainties 
were estimated from the curvature of the surface around the 
minimum, using the standard equation cr - -sJlfW^x^. The lack 
of more than one primary eclipse and the relatively few ASAS 
data points with large errors compared to the NMSU-lm data 
are the main causes of uncertainties on the ephemeris. We have 
no additional constraints on the ephemeris and keep this fixed to 
the above values during the eclipsing binary modelling in Sect. 

s 

With observations of only one complete primary and sec- 
ondary eclipse it is difficult to accurately determine the eccen- 
tricity of the system. However, we estimate this parameter by 
measuring the time of the midpoint of the complete secondary 
eclipse obtained on HJD=2454114 with respect to Tq by mod- 
elling the shape of that eclipse as a 4th order polynomial and 



minimizing the;\f^. We find the timing of the secondary eclipse 
is marginally consistent with a circular orbit for the binary 
within the quoted uncertainties on the ephemeris. Therefore, we 
adopt an eccentricity e - for the system. However, additional 
highly time-sampled observations of multiple primary and sec- 
ondary eclipse minima are necessary to derive a more accurate 
ephemeris, and determine whether there is a small but non-zero 
eccentricity in the system. Figure[T]shows all the small telescope 
photometry phase-folded with the ephemeris defined above. 



2.2. Spectroscopy 

A total of 17 spectra of ASAS J052821H-0338.5 with exposure 
times of 15-30 minutes were obtained with the Fibre-fed Echelle 
Spectrograph (PIES) at the 2.5m Nordic Optical Telescope dur- 
ing January 2-16, 2007. These spectra cover the wavelength 
range 4000-7350 A at a resolution of R ^ 41 000 and have a 
typical signal-to-noise of 30-50. 

The spectra were reduced with the on-line data reduction 
software supplied at the telescope. This package, based on 
PYTHON and PyRAF was especially developed for PIES and 
performs all conventional steps of echelle data reduction, includ- 
ing the subtraction of bias frames, modelling and subtraction of 
scattered light, flat-fielding, order extraction, normalization (in- 
cluding fringe-correction) and wavelength calibration. In addi- 
tion we applied the heliocentric velocity correction to the wave- 
length scale and performed continuum normalization. 

Pour additional 15-minute high resolution (R x 65 000) 
spectra were obtained on April 5, 2007 with the UCLES spec- 
trograph and the SEMPOL polarimeter at the Anglo Australian 
Telescope. These data were reduced with Echelle Spectra 
Reduction: an Interactive Tool (ESPRIT, Donati et al. 11997b . a 
comprehensive reduction package for spectro-polarimetric ob- 
servation. 

The spectrum of ASAS J05282 1-1-0338.5 is a clear example 
of a double-lined spectroscopic binary. Pig. [3] shows a sample 
of the spectrum, obtained near phase 0.64, with well-separated 
components. The pre-main-sequence nature of the object is ev- 
ident from the Ha emission (EW =s 1 A, secondary only) and 
strong Li i 6708 absorption. Except for Ha, no other emission 
lines are seen in the spectrum. Some of the deepest absorption 
lines in the spectrum, such as the Mg b line complex and Na i D 
doublet, show weak emission in the line cores. There is no evi- 
dence of excess continuous emission, such as veiling. 
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3. Analysis 

3.1. Radial Velocity Curve 

To derive the radial velocity curve for the binary system, the 
17 FIES spectra were cross-correlated against the spectrum of 
HD 286264, a single-lined K2IV weak-line T Tauri star, on an 
order-by-order basis using the FXCOR package in IRAF. The 
majority of spectra showed two distinct cross-correlation peaks, 
one for each component of the binary. Thus, both peaks were 
fit independently with a Gaussian profile to measure the veloc- 
ity and velocity errors of the individual components. If the two 
peaks appeared blended, a double Gaussian was fit to the com- 
bined profile. For each of the 17 observations we then deter- 
mined a weighted-average radial velocity for each star from all 
orders without significant contamination by telluric absorption 
features. Here we used as weights the inverse of the variance 
of the radial velocity measurement in each order, as reported 
by FXCOR. The resulting radial velocities were calibrated us- 
ing the velocity offset derived by cross-correlating the template 
spectrum with a radial velocity standard, HD 109358 (Udry et al. 
11999 ). In these data, we find no evidence for a third component, 
since the cross-correlation function showed two distinct peaks. 

To each of the four UCLES spectra, we applied the technique 
of least-squares deconvolution (LSD, see Donati et al. 1997] l. 
This technique, often used in Doppler Imaging, allows us to con- 
struct, with a very high signal-to-noise, a combined-average line 
profile from about 1000 absorption lines in the spectrum. The 
LSD profile obtained in this way shows two distinct absorption 
profiles corresponding to the components of the binary, allowing 
us to measure the individual stellar radial velocities. 

We determined an initial orbit solution to these radial ve- 
locity measurements using the binary code kindly provided by 
D. Gudehus, with the orbital period held fixed at the value de- 
termined from the photometric lightcurve (Sect. [2X2] i. To mini- 
mize the influence of epochs in which the cross-correlation func- 
tions show highly blended peaks, we excluded two epochs that 
occur within 0.05 phase of primary or secondary eclipse. The 
orbit solution was then iterated, and the internal error bars ad- 
justed between iterations, until the final orbit solution yielded a 
reduced of unity. This orbit solution gives a best-fit eccentric- 
ity of e = 0.02 ± 0.02, i.e. formally consistent with a circular 
orbit. We then refined this solution by combining the radial ve- 
locity data with the observed lightcurve (see Sect. |3.3| l. This so- 
lution then yields: a binary mass ratio of Mi /M2 = 1 .03 ± 0.005, 
velocity amplitudes of Ki = 92.4 + 0.33 and K2 = 95.1 + 0.56, 
and a system radial velocity of y = 22.8 ± 0.3 km/s. The final 
radial velocity curve is shown in Fig.|4] 

3.2. Spectroscopic analysis 

In order to study each star individually, we numerically recon- 
structed the spectra of each of the two components with tomo- 
graphic separation. This technique inverts the relation that the 
observed spectrum is a linear superposition of the spectra of the 
individual stars, where the only scaling factor is the luminosity 
ratio of the two components, which depends on their temper- 
atures and radii. Provided that spectra observed at a range of 
phases are available, tomographic separation can efficiently re- 
construct the individual stellar spectra, and is now a standard 
technique for the analysis of double-line spectrosocopic bina- 
ries (see, for example, Bagnuolo & Gies 1991). We applied this 
technique to six spectra, selected on the basis of having the best 
signal-to-noise as well as covering a range of orbital phases. 
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Fig. 4. Radial velocity measurements of the primary star (filled 
circles, dashed line) and secondary star (asteriks, solid line) of 
ASAS J052821+0338.5. A model radial velocity curve for both 
components is overplotted using the ephemeris derived in Sect. 
2.1.2 A binary mass ratio, M1/M2 = 1.03, provides the best fit 



to the data. 



Here, we assumed a luminosity ratio of 2.0 for the two stars, 
which is consistent with the ratio of temperatures and radii ob- 
tained from lightcurve modelling (see Sect. 



3.3 1 



From the individual component spectra we then determined 
the fundamental stellar parameters using the IDL-based spectro- 
scopic analysis package SME (Valenti & Piskunov 1996). SME 
combines radiative transfer calculations with multi-dimensional 
least-squares minimization to determine which parameters (the 
effective temperature Tetf, the gravity log^, the metallicity 
[M/H], the projected rotational velocity vsini, the systemic ra- 
dial velocity v^ad, the microturbulence Vj^i^ and the macrotur- 
bulence Vmac) best describe the observed stellar spectrum. The 
grid of model atmospheres used for radiative transfer are taken 
from Kurucz (1993 ), and atomic line data were obtained from 
the VALD database (Piskunov et al. 1995; Kupka et al. 1999). 
The atomic line data were checked and fine-tuned against a high- 
quality spectrum of the Sun (Kurucz et al. il984 ). 

The two wavelength regions we considered in our spec- 
trscopic analysis are 5880-5920 A and 6000-6200 A. The first 
region contains the Na i D lines, sensitive to both T^ff and 
log g. The second regions contains a wealth of non-blended lines 
of different elements, providing constraints on Teff, \ogg and 
[M/H]. We have no independent leverage on v^ic and v„ac^ and 
therefore assumed v^ic = 1 .6 km/sec (the average value for K- 
type PMS stars found by Padgett 1996) and Vmac = 3.5 km/sec 
(for a T X 5000 K star, see Gray 1992). From this spectro- 
scopic analysis we find for the primary Tgff = 5103 ± 100 K, 
log g -4.0 + 0.1 and [M/H] = -0.2 ± 0.2, and for the secondaiy 
Teff = 4705 + 100 K, log^ = 3.9 ± 0.1 and [M/H] = -0.1 ± 0.2. 
These parameters correspond to spectral types of Kl for the pri- 
mary and K3 for the secondary (Cohen & Kuhi ll979T l. 

The projected rotational velocities of the two components 
were determined by analysing the LSD line profiles (see Sect. 



3. 1 1 of each spectrum. During this process, we took into account 
line broadening from other sources, such as macroturbulence, 
smearing due to the orbital motion of the binary during the ex- 
posure and the instrumental profile. We find an average vsini 
of 24.5 + 0.8 km/sec and 24.5 ± 0.7 km/s for, respectively, the 
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Fig. 5. The above panels show 50 A of the observed spectrum 
of ASAS J05282 1+0338.5 (upper panel), as well as the tomo- 
graphically reconstructed spectra of the primary (middle panel) 
and the secondary component (lower panel), both in their re- 
spective stellar rest frames. The thick solid lines in each panel 
are synthetic spectra calculated with the stellar parameters from 
Table [U 



primary and secondary. These values suggest the stars rotate syn- 
chronously (see Sect.|4]i. 

In addition to the fundamental stellar parameters, we de- 
termined the abundance of lithium from the prominent reso- 
nance Une at 6708 A (EWp,™ ~ 307 mA, EWsec ~ 394 mA). 
After applying corrections of -0.29 (primary) and -0.21 (sec- 
ondary) for NLTE effects (Carlsson et al. 119941) . we obtained 
logn(Li) = {LilH) + 12 = 3.10 + 0.2 for the primary, and 
log«(Li) = 3.35 + 0.2 for the secondary. These values are con- 
sistent with the commonly adopted value of 3.2-3.3 for the pri- 
mordeal lithium abundance. 

Examples of tomographically separated spectra, as well as a 
comparison between the observed and synthetic spectra for each 
star are shown in Fig. |5] A full compilation of the stellar and 
atmospheric parameters is presented in Table [T] 



3.3. Lightcurve modelling 

Together with determining the stellar and orbital parameters 
from spectroscopy we fit the V-band lightcurve and radial ve- 
locity data using the most recent version of the eclipsing binary 
lightcurve modeling algorithm of Wilson & Devinney 



with updates), as implemented in the phoebe code of Prsa & 
Zwitter J2005I I. The code determines the surface gravity and ef- 
fective temperature of each star based on the calculated grav- 
itational potential, assuming full Roche geometry according to 
the formalism of Kopal (e.g. Wilson 1971), and the gravity 
brightening coefficients from Claret ( 2000) . It includes theoreti- 
cal Kurucz model atmospheres to determine intensities over the 
stellar disks. We adopted a linear limb darkening law using coef- 
ficients extrapolated down to 2700 K from Van Hamme (I 1993I I. 
We allowed the code to calculate reflection effects, adopting 
a bolometric albedo of 0.5, typical for a fully convective stel- 
lar envelope. Since the out-of-eclipse variations in the observed 
lightcurve suggest the presence of starspots which are changing 
over time, we include in the modelling only V-band photometry 
obtained with the NMSU-lm and the JGT-lm between January 
9-February 10, 2007. 



To begin, we first applied a traditional rectification procedure 
in which we fit the out-of-eclipse variations with a smooth poly- 
nomial function. The aim of this rectification procedure is to "re- 
move" the spot signal from the lightcurve and to thereby allow 
the lightcurve model to achieve goodness of fit without introduc- 
ing additional spot-fitting parameters that are poorly constrained 
at present. 

Next, in order to maintain control of the lightcurve solution 
and its many free parameters, we performed this fitting in stages. 
In the initial stage we held the orbital parameters fixed at the 
values determined from the double-lined spectroscopic orbit so- 
lution (see Sect. 



3.1 



and fixed the effective temperature of the 
primary component to the value determined from our spectral 
analysis (5103 K; see Sect. 3.2 1. This allowed us to obtain ini- 
tial estimates of the component temperatures and radii, and the 
system inclination. 

With initial values for the temperatures, radii, and inclination 
so determined, we then iteratively improved the solution by first 
allowing the eccentricity and argument of periastron to be fit, 
and then performing a final fit in which all of the orbital and 
component parameters of the system were fit freely (except for 
the ephemeris, see Sect. |2.1.2[ ). The eccentricity we recovered 
for the system, e = 0.007 ± 0.006, is indistinguishable from a 
circular orbit, and we could therefore not constrain the periastron 
angle. Also, if we assume a truely circular orbit (e - 0) we find 
no measurable change in any of the orbital parameters (within 
the limits of the error bars). We have therefore chosen to list 
this system as circular, and report no periastron angle or time 
for periastron passage. The final parameters that result from this 
simultaneous fit to the radial velocities and rectified lightcurve 
are presented in Table [T] 

To explore the influence of, and origin of the out-of-eclipse 
variability, and to better characterize systematic uncertainties 
in the stellar parameters arising from this variability, we have 
also modelled the un-rectified lightcurve by introducing multi- 
ple spots into the fitting analysis. The number of spots, their 
placement on the stellar surfaces, their temperatures and sizes, 
were initially determined manually in order to approximate the 
strong out-of-eclipse variability, then the code was permitted to 
adjust the spot parameters along with all of the other system and 
component parameters. We wish to emphasize that the adopted 
spot solution is neither definitive nor unique. Degenerate starspot 
configurations are consistent with the observed out-of-eclipse 
photometry, and it is not possible to break the degeneracy with 
the current single-band lightcurve. Therefore we present only 
one possible starspot configuration which well mimics the ob- 
served out-of-eclipse variations. Despite the degeneracy of the 
starspot solutions, the lightcurve cannot be adequately fit with- 
out the presence of a large (~ 25° radius) cool spot on the back 
of the primary star We also include a small hot spot near the 
back of the secondary star (i.e. facing away from the primary) 
in order to reproduce the small "bump" just prior to the primary 
eclipse. In Table[T]we summarize the stellar parameters resulting 
from simultaneously fitting the radial velocity and un-rectified 
lightcurve data, for ease of comparison with the parameters re- 
sulting from fitting the rectified lightcurve. The different treat- 
ment of modelling the distribution of starspots did no result in 
any change in the orbital parameters, and these are therefore not 
repeated in Table [T] 

The stellar parameters resulting from the spot modelling pro- 
cedure are overall in good agreement with those resulting from 
the rectification procedure. The formal errors on the stellar pa- 
rameters derived from the above analysis are given in Table [T] 
but these do not include systematic errors inherent to the mod- 
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elling procedures. Still, by comparing the results from the two 
methods (spotted vs. rectified) we find that, despite the obviously 
strong eff'ect of spots on the lightcurve, a detailed treatment of 
the spots is evidently not critical for most of the stellar parame- 
ters of interest. The notable exceptions are the temperature and 
radius of the secondary. These parameters depend on the eclipse 
depths, which are affected by the manner in which the out-of- 
eclipse variations are modelled. 

The temperature of the secondary is determined relative to 
the temperature of the primary (which we hold fixed at 5103 K) 
by virtue of the temperature ratio, which is principally deter- 
mined from the ratio of eclipse depths. Between the two methods 
of analysis, we find a difference of ~ 150 K in the inferred tem- 
perature of the secondary, a difference that is significantly larger 
than the formal uncertainty of ~ 25 K in the lightcurve model 
(Table[T]l. The temperature resulting from the rectified analysis is 
more consistent with the value from spectroscopic analysis of the 
tomographically separated spectra (Sect. 3.2 1, thus in what fol- 



lows we adopt this solution as the preferred one. This lightcurve 
model is overplotted on the phase-folded data in Fig. [6] 

The radii obtained from the two light-curve solutions also 
show larger differences than expected from the formal errors. 
Even so, the sum of the radii is consistent in both solutions to 
within ~ 2%. Due to the circular orbit, the primary and sec- 
ondary eclipses are of equal durations and thus it is not possible 
to determine the individual radii independently from a single- 
band light curve. While the sum of the radii is very well con- 
strained by the total duration of the two eclipses, almost any 
combination of component radii that satisfy the sum will yield a 
good fit to the light curve. Of course, the diff'erent component 
radii will in general result in very different flux ratios (since 
the temperature ratio is well constrained). Thus, we caution that 
additional observational constraints (e.g. flux ratios at multi- 
ple wavelengths) will be needed to firmly establish the individ- 
ual radii accurately. The spectroscopically derived temperatures 
cannot by themselves constrain these flux ratios, since tomo- 
graphic separation requires us to assume a luminosity ratio (see 



Sect. 3.2 1. We have therefore iteratively worked towards a solu- 
tion for the temperatures and radii that is consistent with both 
the spectroscopic and the photometric models. 



Table 1. Properties of ASAS J052821 -1-0338.5. For each com- 
ponent we list the atmospheric parameters: effective tempera- 
ture (Tgff), surface gravity (logg), rotational velocity (vsin/), 
metallicity {[M/H]) and lithium abundance (logn(Li)). We also 
present the systems ephemeris: the ephemeris reference epoch 
(To) period (P), eccentricity (e), systemic radial velocity (y), in- 
clination (/) and projected semi-major axis (a sin/)- Since the 
system's eccentricity is indistinguishable from truly circular, we 
choose e = and list no periastron angle or time for periastron 
passage. Finally, we list the parameters of our solutions from 
lightcurve modelling with and without spots: the velocity ampli- 
tude (K), stellar mass (M), radius (R), luminosity (L) and effec- 
tive temperature (Teff) 



Parameter 



Primary 



Secondary Unit 



Stellar parameters from spectral analysis 



V sin I 

[M/H] 

logn(Li) 



5103 4 

4.0: 

24.5 : 

-0.2: 
3.10: 



100 
0.1 
0.8 
0.2 
0.2 



4705 4 

3.9: 
24.5: 
-0.1 : 

3.35: 



100 
0.1 

:0.7 
:0.2 
:0.2 



K 



km s 



Parameters from rectified lightcurx'e modelling 



To 


HJD 2454110.8500 ±0.005 




P 


3.8729 ± 


0.0002 




d 


e 











7 


22.8 4 


0.3 




km s 


i 


83.7 4 


1.0 




deg 


a sin i 


0.0668 ± 


0.0002 




AU 


K 


92.41 ±0.33 


95.08 


±0.56 


km s 


M 


1.375 ±0.011 


1.329 4 


0.008 


Mq 


R 


1.83 ±0.01 


1.73 


±0.01 


Rq 


L 


2.05 ±0.16 


1.38 


±0.11 


Lq 


Tcff 


5103* 


4751 ±26 


K 


^ogg 


4.05 ± 0.01 


4.08 


±0.01 





Parameters from EB modelling with spots 



M 

R 

L 

log 5 



1.387 ±0.017 
1.84 ±0.01 
2.06 ±0.16 
5103* 
4.05 ± 0.02 



1.331 ±0.011 
1.78 ±0.01 
1.28 ±0.10 
4599 ± 23 
4.06 ±0.01 



Mq 
Rq 
Lq 
K 



3.4. Distance, physical iocation and kinematics 

We estimated the distance to ASAS J052821+0338.5 by com- 
paring distance moduli from the available broad-band magni- 
tudes. For this analysis, we have assumed relatively large uncer- 
tainties on the observed magnitudes (0.15'"), because we have 
no possibility of knowing whether these are truly out-of-eclipse 
values. Using the components' surface temperatures of the pre- 
ferred solution from lightcurve modelling (5103 K and 4751 K) 
and the standard color indices tabulated by Kenyon & Hartmann 
( |1995t . we find that the distance moduli of the available col- 
ors converge on a distance of 280 + 30 pc and an extinction of 
Ay - 0.6 + 0.3™. The distance is mostly constrained by the 
(almost) extinction-insensitive J, H and K magnitudes, while 
the extinction depends mostly on the B and V magnitudes. The 
quoted error margins are dominated by the uncertainties of the 
observed magnitudes. 

In addition to the above, we compared the observed broad- 
band magnitudes of ASAS J052821 +0338.5 wifli model stellar 
atmosphere spectra from Castelli & Kurucz (2003), using the 
SED fitting routines described in Robitaille et al. (2007; modi- 
fied to incorporate stellar atmosphere models). Assuming a tem- 
perature of 5000 K, which is the luminosity-weighted average 



* Assumed parameter for EB modelling 



temperature of the binary, we find the broad-band magnitudes to 
be compatible with an extinction of Ay x 0.65™. We also mea- 
sured the equivalent width of the narrow interstellar absorption 
components in the Na i Dl line. We find two absorption compo- 
nents with a combined EW of about 200 mA, which translates 
to Av ~ 0.25 ± 0.1™ (Munari & Zwitter fTWt . The interstel- 
lar Na I D absorption is partly filled-in by telluric emission, and 
therefore this value should be considered a lower limit to Ay. We 
notice that the strongest absorption component (with an EW of 
~ 150 mA) has a central velocity equal to the systemic veloc- 
ity of the system, suggesting that a large part of the extinction 
originates from regions close to the binary. 

Following the definitions of Warren & Hesser (1977]l, 
ASAS J052821 +0338.5 is located in the direction of the widely 
spread Orion OB la subassociation. The system is about 3^ de- 
grees away from the much more compact Orion OB lb subas- 
sociation and about 2 degrees from the 25 Ori subassociation, 
a compact and kinematically distinct subgroup of Orion OB la 
recently identified by Bricefio et al. (|2007 ). Orion OB la is lo- 



H. C. Stempels et al.: The Pre-Main-Sequence Eclipsing Binary ASAS J052821+0338.5 7 




Fig. 6. The top panel shows the phase-folded and rectified diff'erential V-band lightcurve obtained with the NMSU-lm and JGT-lm 
telescopes. The lightcurve has a precision of 0.005-0.01™ and a sampling of < 1 min. The data are phase-folded with a period, 
P - 3.8729 days, and epoch, Tq - HJD 2454110.8500, as well as rectified by fitting a smooth polynomial to the out-of-eclipse 
variations. The smooth grey curve shows the calculated binary lightcurve that has the best agreement with the data. The residuals of 
this curve with the data are shown in the bottom panel. 



cated at a distance of 330 + 15 pc, while Orion OB lb is further 
away, at a distance 440 ± 20 pc (Brown et al. .l994, de Zeeuw et 
al. 119991 Hernandez et al. 120051) . There is little diff'erence in the 
kinematic distribution of Orion OB la and Orion OB lb. The av- 
erage radial velocities of Orion OB la and Orion OB lb, 23.8±0.7 
km/sec and 23.1 + 1.4 km/sec respectively (Morrell & Levato 
I1991I I. and also their tangential velocity distributions are similar 
= 0.44 + 0.25 mas/yr, yU^ = -0.65 ± 0.25 mas/yr. Brown et al. 
119981 de Zeeuw et al. 1999). Orion OB la is commonly assumed 
to be 11.4 + 1.9 Myr old, signific antly o lder than the 1.7 ± 1.1 
Myr old Orion OB lb (Brown et 



The systemic radial velocity we determined for 
ASAS J052821H-0338.5 is 22.8 + 0.3 km/sec, and the tan- 
gential velocity from the NOMAD catalog is - -1.6 + 2.8 
mas/yr, and - -4.0 + 4.0 mas/yr These kinematics combined 
with the estimated distance and physical location suggest that 
ASAS J052821-h0338.5 is most probably a member of the 
~ 1 1 Myr old Orion OB la subassociation. 



3.5. Evolutionary status 

Several characteristics of the two binary components, including 
their large radii, their cool temperatures, and their primordial 
lithium abundances, as well as the proximity to the Orion OB la 
subassociation strongly indicate that ASAS J052821+0338.5 is 
a young, pre-main sequence eclipsing binary. 

In Figure |7] the structural properties of the binary are shown 
in relation to other PMS and main-sequence eclipsing bina- 
ries, as well as theoretical stellar evolution models from BarafFe 
et al. ( 119981 1. Both components have larger radii than main- 
sequence objects of the same mass, as expected for stars that 
are still evolving onto the main sequence. However, their large 
observed radii do not alone imply a pre-main-sequence evolu- 
tionary state. A typical, post-main-sequence, 1.4M0 star with 
solar metallicity reaches a radius similar to that of the primary 
component of ASAS J052821-H0338.5 at an age of ~ 2 Gyr. 
Yet, at about 6000 K, such a star is considerably hotter than 
either component of ASAS J052821-H0338.5. A 1.4Mo post- 
main-sequence star will cool to a temperature comparable with 
ASAS J052821H-0338.5 (~ 5000 K) as it evolves through the 
red giant phase, but this necessarily coincides with a further in- 
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Fig. 7. Mass versus radius for known main sequence (open cir- 
cles) and pre-main sequence (grey solid circles) eclipsing bina- 
ries. The components of ASAS J05282 1-1-0338.5 are shown as 
black solid circles. Isochrones from Baraffe et al. d 19981 1 with 
ages of 3 Myr, 6 Myr, 10 Myr, 13 Myr, and 300 Myr are over- 
plotted. 



crease in radius to a value much larger than what is observed in 
the binary. Thus, the cool temperatures of the binary components 
combined with their ~ 1 .4 M0 masses and < 2R0 radii suggest 
ASAS J052821 +0338.5 is a pre-main sequence object. 

This is exemplified in the modified Hertzprung-Russel (HR) 
diagram shown in Fig. [8] In this figure, stellar radius is a proxy 
for luminosity, as it is a direct measurement from the eclips- 
ing binary analysis, and it relates to luminosity through the 
temperature. Evolutionary tracks and isochrones for the Baraffe 
et al. (1 19981 1 models are shown in addition to the position of 
ASAS J052821-I-0338.5 and the other known PMS eclipsing bi- 
naries. We interpolated the evolutionary tracks at the masses of 
the primary (black solid line) and secondary (black dashed line) 
components. ASAS J052821+0338.5 is clearly positioned in the 
PMS phase of evolution in this parameter space when compared 
to empirical data and to this set of models. 

A further indication of the PMS evolutionary state of 
ASAS J052821 -1-0338.5 is the detection of primordial Uthium 
in both components. The rate of lithium depletion in late-type 
stars depends on the depth of the convection zone (and thus, 
on temperature), and the primary is therefore expected to have 
the slowest lithium depletion rate. The primary, with an effec- 
tive temperature of 5100K, will deplete lithium well below the 
primordial level after about 150-200 Myr (Sestito & Randich 
12005 In conclusion, only very young stars can have bloated 
radii, cool temperatures, and primordial lithium. 

Finally, although determining an absolute age of 
ASAS J052821 -1-0338.5 would involve a comparison with 
largely uncalibrated stellar evolutionary models, we do note 
that the fundamental parameters of ASAS J052821 -1-0338.5 are 
consistent with the 10-13 Myr solar-metallicity isochrones of 
Baraffe et al. (1998). This is in agreement with the age of the 
Orion OBla subassociation (1 1.4 + 1.9 Myr, see Sect. 3.4 1. 



■D 
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4. Discussion 

The photometric and spectroscopic observations we obtained of 
ASAS J052821-1-0338.5 identify this system as a new pre-main- 



0.5 - 



Fig. 8. The modified Hertzprung-Russel diagram showing ra- 
dius versus effective temperature for the components of 
ASAS J052821 +0338.5 (black soHd circles) and other PMS EBs 
in Orion (grey solid circles). Evolutionary tracks (1.4, 1.3, 1.2, 
1.1, 0.7 M0) and isochrones (3, 6, 10, 13, 16, 300 Myr) from 
Baraffe et al. 1998 are overplotted as grey lines. The black solid 
line shows the evolutionary track interpolated at the mass of 
the primary star (1.375 Mq). The black dashed line shows the 
evolutionary track interpolated at the mass of the secondary star 
(1.329 Mo). 



sequence eclipsing binary with a period of 3.87 days. The masses 
and radii of the components place the system in an as yet unsam- 
pled region of the mass-radius and HR diagrams, and thereby 
provide a new pair of reference points for PMS stellar evolution- 
ary models. Both components of the system are more massive 
than all other members of known PMS eclipsing systems, except 
for TY CrA B, the 1.64 M© and ~ 3 Myr old companion to a 
main-sequence B-type primary (Casey et al. ll998l l. 

The physical location and kinematic properties of 
ASAS J052821+0338.5 suggest this system to be a mem- 
ber of the ~ 11.4 ± 1.9 Myr Orion OBla subassociation. This 
age is also compatible with the position of the two components 
in the mass-radius diagram shown in Fig. [7] where both stars fall 
close to the 13 Myr isochrone of the Baraffe et al. ( 1998) evolu- 
tionary models. At this age, the two 1.35 Mq stars have just left 
the fully-convective Hayashi stage. The temperature of the two 
components will then increase at almost constant luminosity, 
first slightly growing in radius before their final contraction 
towards the main sequence, settling at a final temperature of 
~ 6250 K (for a more detailed description of pre-main-sequence 
evolution, see for example Palla & Stabler 119931 ). We also 
find a that the locations of the stars in the HR diagram agree 
quite well with the predicted 13 Myr evolutionary tracks of 
Baraffe et al. (1998, see Fig. [8]). Although there appear to be 
small differences between the observed radii and the predicted 
values, this does not indicate any inconsistency since the only 
well-constrained parameter is the sum of the radii (see Sect. 
3.3 1, and full consistency of both stars with the predicted tracks 
requires only small antagonal adjustements (of about 3%) to the 
stellar radii. Therefore, a critical next step to further improve 
these parameters will be a multi-band lightcurve analysis with a 
more detailed spot treatment. 

Most of the known PMS eclipsing binaries are located in 
Orion. ASAS J052821+0338.5 and RXJ 0529.4+0041A are 
thought to be located in the Orion OBla subassociation (Covino 
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et al. I2OOOI Covino et al. |2004] |. VI 174 Ori is consid ered 
to be a member of Orion OBlc (Stassun et al. 2004) and 
2MASS J05352 184-201 30546085, JW 380 and Par 1802 are 
found in the Orion Nebula Cluster (Stassun et al. 2006, Irwin 
et al. I2007b1 Cargile et al. 2007 ). These OB subassociations in 
Orion form an evolutionary sequence, with ages of ~ 1-12 Myr 
This age range covers a key period in star formation when stars 
experience rapid evolution in both radius and temperature. The 
components of the PMS eclipsing binaries in Orion have masses 
between 0.03 and 1.35 M0 and have well-determined radii and 
temperatures. As such they form an easily accessible and ho- 
mogenous sample of late-type PMS stars that can be used to per- 
form relative calibrations of PMS stellar evolutionary models. 

As described in Sect. |2.1.1| we detected considerable 
(~ 0.15™) out-of-eclipse variations, occuring around the same 
orbital phase throughout our four-month campaign of photo- 
metric observations, not unlike the 'photometric wave' that is 
characteristic of the chromospherically active RS CVn binaries 
(see the review by Rodono I1992l l. Like in an RS CVn binary, 
the components of ASAS J052821-I-0338.5 have relatively short 
rotational periods and extensive convective envelopes, which 
help to drive a stellar dynamo and enhance magnetic activity. 
The photometric variations may then be explained by a large 
cool spot (or region of spots) located on the back of the pri- 
mary, rotating in and out of view in phase with the binary pe- 
riod. However, this configuration requires the primary to rotate 
synchronously with the binary orbit. According to the theory 
of Zahn & Bouchet (1989), tidal breaking during the Hayashi 
phase, when stars are fully convective, is a very efficient mech- 
anism for inducing orbital synchronization and circularization 
in short-period late-type binaries. Once the convective envelope 
starts decreasing, tidal breaking loses its efficiency, and further 
stellar contraction may lead to (temporary) departures from syn- 
chronous rotation. ASAS J052821+0338.5 seems to follow this 
pattern. The orbit is almost perfectly circularized, and the ob- 
served rotational velocities of the primary and secondary are 
24.5 ± 0.8 km/sec and 24.5 ± 0.7 km/s, respectively. Since the 
expected synchronous rotational velocities are 24.0 km/sec and 
23.4 km/sec, the primary is most likely in synchronous rotation 
with the binary orbit, while the smaller secondary may rotate 
slightly super-synchronously. 

While an extended cool spot on the back of the primary 
can fully explain the observed out-of-eclipse photometric vari- 
ations, we did not detect any signatures of a cool spot in our 
current high-resolution spectra, such as changes in the observed 
line-ratios or the appearance of molecular features. However, the 
brightness contrast between the spotted and the unspotted areas 
is large, and since the signal-to-noise of our spectra is moderate 
it is unclear whether such an effect would be detectable. An alter- 
native explanation to the out-of-eclipse variability might be peri- 
odic eclipses of the system by a warped circumbinary disk, or by 
circumstellar material in co-rotation with the binary. Such con- 
centrations of circumstellar matter localized around the stable 
co-rotating Lagrangian points (L2 and L3) have been detected in 
the close PMS binary V4046 Sgr (Stempels & Gahm l2004l) . We 
are currently pursuing multi-band photometry of the system in 
order to determine the nature of the out-of-eclipse variations, as 
well as to further improve the orbital and stellar parameters of 
ASAS J052821H-0338.5, in particular the temperatures and radii 
of the components. 
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